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Today  we  would  like  to  discuss  some  problems  involving  Long-range  reverberation 
in  shallow  water. 


Fig.  1 


The  title  of  the  talk  is  “Nonreciprocity  of  Long-range  Reverberation  in  Wedged 
Continental  Shelf’.  The  talk  is  divided  into  four  parts.  First  we  would  like  to  talk 
a  little  bit  about  the  background  of  current  research  on  long-range  reverberation  in 
shallow  water.  The  problem  seems  to  fall  into  a  “dilemma”  (endless  loop)  with  no  way 
out.  Second,  we  introduced  an  averaged  angular  power  spectrum  method  for  long- 
range  reverberation  in  the  Pekeris  model  to  show  how  we  could  bypeiss  such  ‘dilemma’. 
Based  on  a  resulted  transformation  relation  between  arbitrary  angle  dependence  of 
bottom  scattering  and  the  range  dependence  of  reverberation  in  shallow  water,  some 
experimental  data  on  bottom  scattering  at  small  grazing  angles  and  low  frequencies 
will  be  given  in  the  third  part.  In  the  last  part,  using  the  WKB  approximation  to 
the  adiabatic  normal  mode  theory,  we  extend  the  problem  into  wedged  homogeneous 
shallow  water.  The  result  shows  that  the  monostatic  reverberation  intensities,  obtained 
at  two  terminals  with  different  depths,  would  not  be  reciprocal.  Because  of  time 
limitations,  the  results  are  presented,  rather  than  developed. 
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I.  The  status  of  current  research  on  long-range  reverberation  in  shallow 
water 


Fig.2 


Generally  speaking,  the  ‘long-range’  means  low-frequency,  and  the  range  is  more 
than  several  hundred  times  of  water  depth.  Here  are  two  typical  long-range  rever¬ 
beration  decay  cmves  for  two  frequencies(1000Hz,  1600Hz;  1/3  OCT),  obtmned  from 
a  shallow  water  with  the  depth  of  29m  using  a  explosive  source.  The  ordinate  is  the 
relative  reverberation  level;  The  abscissa  is  the  range-depth  ratio.  If  a  water  depth  is 
100m,  the  ratio  value  of  400  corresponds  that  the  reverberation  signal  comes  from  40 
km  far  away..  These  reverberation  curves(see  Fig.l)  were  obtained  in  a  Pekeris  shallow 
water  ,  i.e.,  the  simplest  shallow  water  model. 


Fig.3 


In  the  Pekeris  model,  the  effect  grazing  angle  of  sound  propagation  can  be  ex¬ 
pressed  as  Eq.(l),  where  Q  is  the  bottom  reflection-loss  peirameter  at  small  grazing 
angle,  defined  as 


-in\v{9)\  =  qe 
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Here  axe  some  Zhou’s  experimental  results(from  vertical  coherence  meeisurements  of 
sotmd  propagation),  compared  with  Zhou  and  Smith  theory.  The  results  show  that  at 
long-range,  the  grazing  angles  of  sound  interacting  with  bottom  are  very  small.  For 
example,  at  a  distance  of  400  times  of  water,  the  main  energy  come  from  angles  that 
are  less  2.4°  for  800Hz,  less  than  2°  for  1600  Hz.  Due  to  mode  stripping,  the  farther 
the  distance,  the  smaller  the  grazing  angles  of  the  effective  sound  propagation.  What 
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about  the  bottom  scattering  at  such  small  grazing  angles  and  low  frequencies? 


Fig.4 


The  upper  curves  of  Fig.4  is  a  copied  from  Urick’s  book.  At  small  grazing  angles  , 
no  data.  The  follovring  is  quoted  from  the  proceedings  of  1983  Shallow  Water  Acoustics 
Workshop: 

“Grazing  angles  of  primary  importance  to  shallow  water  applications  range  from 

about  30°  to  near  0°  with  the  smaller  angles  being  more  important .  To 

date,  NOVOCEANO  has  not  been  able  to  report  either  bottom  reflectivity  or  bottom 
backscattering  values  for  any  shallow  water  area.” 

Nine  years  have  past,  as  we  know,  still  no  any  progress  in  this  area.  The  multipath 
in  shallow  water  and  the  unclear  mechanisms  of  bottom  scattering  at  low  frequency 
make  the  computation  of  long-range  reverberation  very  complex.  The  theoretical 
acousticians  have  no  adequate  and  reliable  data  base  for  bottom  scattering  to  develop 
a  practical  theoretical  model.  Experimental  Mousticians  say  we  have  no  suitable  theo¬ 
retical  reverberation  formula  to  use  for  measuring  bottom  scattering  strength  at  small 
grazing  angle.  (  Because  the  signiflcant  diiflculty  is  that  if  the  assumption  used  in  ex¬ 
tracting  the  scattering  coeflicients  from  the  received  reverberation  has  no  enough  data 
as  a  base,  the  resulting  coefficients  will  not  represent  valid  environmental  parameters 
which  should  be  independent  of  how  the  data  were  obtained  and  analyzed.)  Thus  the 
problem  of  long-range  reverberation  in  shallow  water  seems  to  fall  into  an  endless  loop 
with  no  way  out. 

Some  groups  eire  working  on  developing  theoretical  model  of  bottom  scattering 
at  low  frequency,  and  trying  to  get  its  angular  and  frequency  dependence.  A  pressing 
problem  is  that  how  to  judge  the  validity  of  their  theoretical  results  at  low  frequency 
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and  low  grazing  angle.  (  Currently  available  data  are  limited  to  high  frequencies  or 
larger  grazing  angles.) 

Next  we  would  like  briefly  to  introduce  a  averaged  angular  power  spectrvun  method 
for  shallow  water  long-range  reverberation,  and  to  give  some  data  on  bottom  scattering 
at  small  grazing  angles  and  low  frequencies. 

II.  An  averaged  angular  power  spectrum  method  for  long-range  reverber¬ 
ation 


Fig.  5 


For  a  point  source  in  shallow  water,  applying  the  WKB  approximation  to  the 
normal  mode  expression(Brekhovskih,Zhou),  or  using  the  ray  method(Smith)  or  the 
flux  method(Weston)  the  averaged  sound  field  intensity  can  be  expressed  as  Eq.(2) 
where  0(2o)^and  6{z)  are  the  equivalent  ray  angle  of  normal  mode  at  the  source 
depth  or  receiver  depth,  S  is  mode  cycle  distance.  ln\V{9)\  is  the  bottom  reflection- 
loss. 

We  call  laps  the  angular  power  spectrum.  Except  the  cylindrical  spreading  and  the 
medium  absorption,  the  sound  velocity  profile  and  the  boimdary  condition  together 
compose  a  angular  filter  of  shallow  water  sound  propagation.  The  angular  expression  of 
sound  propagation  can  naturally  be  connected  with  the  classical  expression  of  bottom 
scattering  in  the  angular  domain.  For  a  given  signal  duration,  there  are  a  lot  of 
scattering  signal  simultaneously  return  to  a  receiver.  The  bottom  scattering  looks  like 
a  stochastic  filter.  The  reverberation  can  be  treated  as  a  angle-weighted  process.  The 
net  analog  of  shallow  water  reverberation  can  be  expressed  as  Eq.  (3).  The  forward 
transmission  net,  bottom  scattering  net  plus  backward  transmission  net. 

Where  the  M{9,<f>)  is  the  bottom  scattering  coefficient  for  plane  wave.  9  is  the 
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grazing  angle  of  incident  mode-ray,  and  <{>  is  the  scattering  angle.  .A  is  a  scattering  area. 
The  validity  of  a  theoretical  result  from  this  expression  will  depend  on  the  validity  of 
scattering  coefficient  M{d,<f>).  Unfortunately  as  above-mentioned,  at  low  frequencies 
and  low  grazing  amgles  there  is  no  data  or  general  expression  about  the  It 

would  be  very  complex  function  of  bottom  surface  roughness,  sediment  type,  sediment 
inhomogeneity,  angle,  frequency  and  so  on.  It  still  seems  no  way  out. 

Here  we  use  a  little  trick  to  bypass  the  “endless  loop” . 


Fig.6 


Here  is  an  arbitrary  angle  dependence  of  bottom  scattering,  ignoring  the  concrete 
analytical  expression,  i.e.,  no  matter  what  is  the  main  mechanism  of  bottom  scatter¬ 
ing,  surface  roughness,  volume  inhomogeneity  or  whatever.  General  speaiking,  The 
scattering  coefficient  M{d,  <t>)  is  a  slowly  changed  function.  The  smaller  the  angle  9  or 
<l>,  the  smaller  the  scattering  strength.  If^we  consider  backward  scattering  is  reciprocal, 
i.e.,  M{9,4>)  =  M{(i),9),  then  at  small  grazing  ^mgles  any  scattering  function  can 
always  be  expressed  2is  Eqs.(4-5). 

We  introduce  this  phenomenological  expression  only  as  a  “bridge”  between  long- 
range  reverberation  and  bottom  scattering,  ignoring  the  concrete  analytical  expression. 
Next  we  will  see  that  it  is  this  little  trick  to  help  us  to  bypass  the  above-mentioned 
“endless  loop”. 


Fig.7 


In  the  Pekeris  model  the  angular  power  spectrum  is  expressed  by  Eq.(6). 

From  Eos.  (3-6),  we  get  an  expression  for  long-range  reverberation  as  Eqs.(7-8). 
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(The  reverberation  time  tj  =  ^,  c-sound  velocity.) 

The  result  shows  that  the  characteristics  of  long-range  reverberation  between 
distance  r,-i  and  mainly  depend  on  the  bottom  scattering  indices  Uj  in  the 
region  between  Qi^x  and  9i.  The  transformation  between  the  range  dependence  of 
reverberation  level  RL  and  the  angle  dependence  of  bottom  backscattering  strength  is 
shown  in  Fig.  7. 

From  this  transformation  relation  the  scattering  indices  /Xj  and  rij  can  be  derived 
from  two  data  on  the  experimental  curve  of  reverberation  level  near  rj.  Then  we  obtain 
the  bottom  backscattering  strength  at  the  angle  9i  (see  Eq.(l)). 

III.  The  reverberation-derived  bottom  scattering  strength  at  small  grazing 
angels. 


Fig.  8 


The  bottom  scattering  strength  for  frequency  band  of  800Hz-AkHz  md  grazing 
angles  of  2°-8°  ,  shown  here,  is  derived  from  at-sea  experimental  reverberation  data. 

For  the  reliability  of  results  only  those  experimental  long-range  reverberation  data 
were  used  for  extracting  bottom  scattering  that  were  about  9  dB  higher  than  the 
environmental  noise.  (For  comparison,  the  bottom  scattering  at  larger  angles  obtained 
by  other  researchers  from  deep  water  measurements  is  included  with  dote  or  dashed 
lines.) 


Fig.  9 


Here(Fig.  9)  is  the  bottom  scattering  strength  derived  from  another  reverberation 
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experiment. 


Fig.  10 


Russian  acousticians  Ivakin  and  Lysanov  compared  their  theoretical  bottom  scat¬ 
tering  model  with  our  experimental  data(Fig.  10),  the  discrepancy  does  not  exceed 
3  dB  at  any  of  the  investigated  angles  and  frequencies.  Our  results  have  two  specific 
features  which  are  different  from  high  frequency  results:  bottom  scattering  has 
more  strong  frequency  dependence,  and  decrease  much  more  rapidly  as 
the  grazing  angle  is  decreased.  At  very  small  angles,  the  bottom  scattering  index 
could  be  larger  than  Lambert’s  law  scattering  i.e.,  n  >  2. 

rV.  Long-range  reverberation  in  wedged  homogeneous  continental  shelf 


Fig.  11 


Next  we  extend  the  results  of  long-range  reverberation  to  shallow  water  of  variable 
depth,  (we  consider  the  case  of  propagation  directly  up  or  down  slop  from  the  source 
) 

If  the  water  depth  and  sound  velocity  profile  vary  extremely  slowly  with  hori¬ 
zontal  distance,  by  neglecting  any  coupling  between  normal  modes,  Pierce  obtained 
a  adiabatic  normal  mode  expression  for  sound  propagation  as  Eq.(9)  (a  cylindrical 
coordinate  has  been  assumed): 

where . 

The  adiabatic  mode  theory  means  that  the  local  modes  adapt  to  local  environ¬ 
ment.  So  the  grazing  angle  of  equivalent  ray  of  normal  modes  satisfies  Eqs. (10-11). 

For  a  wedged  homogeneous  continental  shelf  the  water  depth  can  be  expressed  by 
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Eq.(  12). 

Applying  WKB  approximation  to  the  Pierce’s  expression,  the  averaged  sound 
intensity  in  homogeneous  wedged  shallow  water  can  be  expressed  as  Eq.(13)  . 

The  averaged  angular  power  spectrum  in  wedged  water  is  expressed  by  Eq.(14). 


Fig.  12 


Using  above-mentioned  angular  power  spectrum  method,  we  get  a  averaged  sound 
intensity  in  the  most  interested  three-half  law  field  area  as  Eq.(15).  The  averaged 
reverberation  intensity  is  expressed  by  Eq.(16). 

The  sound  propagation  is  reciprocal  if  we  exchange  the  positions  of  the  source  and 
receiver.  But  the  monostatic  reverberation  intensities,  obtained  at  two  terminations 
with  a  depth  of  H\  or  H2,  would  not  be  reciprocal  (see  Eq.(17)). 

Here  N  is  the  angular  index  of  bottom  scattering  defined  as  Eq.(18).  If  there 

are  two  identical  active  sonars,  but  the  water  depth  of  their  locations  is  different,  the 

♦ 

reverberation  interference  to  the  sonar  at  shallower  area  would  be  much  smaller.  For 
example,  if  Hi  =  50m,  H2  =  200m,  N  ■=  2  (so-call  the  Lambert’s  law  scattering), 
then  the  reverberation  interference  at  shallower  area  would  be  12  dB  less. 

At  last  year’s  Shallow  Water  Acoustic  Workshop  held  at  Woods  Hole,  a  report 
firom  NUSC,  New  London  mentioned  that  the  reverberation  in  continental  shelf  is  not 
reciprocal.  The  authors  of  this  talk  hope  their  theoretical  result  could  be,  in  some 
degree,  helpful  to  explain  the  NUSC  at-sea  experiments. 

Thanks  ! 
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Nonreciprocity  of  Long-range  Reverberation 
in  a  Wedged  Continental  Shelf 

J.  Z.  Zhou,  X.  Z.  Zhang,  P.  H.  Rogers  and  D.  H.  Guan 


1.  Background  (The  shallow  water  reverberation  problem  seems 
to  fall  into  “an  endless  loop’’  with  no  way  out.) 

2.  Averaged  angular  power  spectrum  method  for  long-range 
reverberation  in  shallow  water. 

3.  Transformation  relation  — >  Bottom  scattering  at  small  grcizing 
angles,  derived  from  experimental  reverberation  data. 

4.  The  WKB  approximation  — >  adiabatic  mode  theory. 
Reverberation  in  a  wedged  shallow  water  is  not  reciprocal. 
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-in\v(e)\  =  Qe 

Q  -  the  bottom  reflection  loss  parameter  at  small  grazing  angle  9 
The  Vertical  coherence  of  sound  propagation  in  Pekeris  model(Smith,  Zhou) 
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FIG.  8-22.  Low-frequency  deep-water  bottom  back- 
scattering  strength  as  reported  in  four  sets  of 
measurements. 


Shallow  Water  Acoustics  Workshops  Proceedings  (1983)  : 

“Grazing  angles  of  primary  importance  to  shallow  water  applications  range 

from  about  30°  to  near  0°  with  the  smaller  angles  being  more  important . 

To  date,  NOVOCEANO  has  not  been  able  to  report  either  bottom  reflectivity 
or  bottom  backscattering  values  for  any  shallow  water  area.  ” 
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Brekhovskih,  Zhou 
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hpsi^i  r,  z;  zq)  -  the  angular  power  spectrum; 

Boudary  condition  +  sound  velocity  profile  — >  a  angular  filter  of  sound 
propagation. 
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where  the  AI{0.  4>)  is  the  scattering  coefficient  of  plane  wave  from  bottom. 

M{9,  (i>)  =  ? 

It  must  be  a  very  complex  function  of  bottom  surface  roughness,  sediment  type, 
sediment  inhomogeneity,  angle,  frequency  and  so  on.  Unknown  !  No  way  out  ? 
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Arbitrary  angle  dependence  of  equivalent  bottom  scattering: 


Reasonable  assumptions  at  small  grazing  angles 

1)  Decreasing  function  of  angle  with  decreasing  gra,zing  angle; 

2)  Reciprocal  when  incident  angle  and  backward  scattering  angle  exchange, 
i.e.,  M(9,4>)^M(<i>,e) 


where 


ifxG 

'  \  0,  otherwise. 


(5) 


RL(dB) 


In  the  Pekeris  model  the  angular  power  spectrum 


IaA^,r)  =  e^^ 

From  Eqs.(3-6),  we  get  an  expression  for  long-range  reverberation  as 

N 
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IVansformation  relation  between  the  range  dependence  of  RL 
and  the  angular  dependence  of  bottom  scattering: 
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FIG.  1.  Intensity  of  sound  scattering 
by  tic  bottom  vs  grazing  angle  at  vari 
ous  frequencies:  1)  0,8  kHz;  2)  1.0; 
3)  1.6;  4)  2.5;  5)  4  kHz. 
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FIG,  2,  Intensity  of  sound  scattering 
by  the  bottom  vs  frequency  at  various 
grazing  angles:  1)  2.5*;  2)  5*. 


In  shallow  water  of  variable  depth,  the  an  adiabatic  normal  mode  expres¬ 
sion  (Pierce)  : 
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The  adiabatic  mode  theory  means  that  the  local  modes  adapt  to  local 
environment.  So  the  grazing  angle  of  equivalent  ray  of  normal  modes  satisfy 

km{r)  =  k{r)cos{em{r))  (10) 


ff(r)sin^„,(r)  =  .ff(O)sin0,„(O)  (11) 

For  a  wedged  homogeneious  continental  shelf  the  water  depth  can  be  ex¬ 
pressed  by 

H{r)=H{0){l  +  ar)  (12) 

Appling  WKB  approximation  we  have  the  averaged  sound  intensity  : 


/(r)  = 


1^^  [L  - 


I  Vi^^n)  I 


H{0)r  J  [Jo  H{r) 

The  averaged  angular  power  spetrum 


tan(9^(r))dr 


d9m{r)  (13) 


dr 


#  11 


(14) 


The  averaged  intensity  of  sound  propagation  is: 


I{r)  = 


27r 


Q[/f(0)  +  fl'(r)]  r3/2 

The  averaged  reverberation  intensity  is: 


(15) 


Q  ^■(0)+g(^)]fl+n/2^  r^2+n/2) 


(16) 


Sound  propagation  is  reciprocal  if  we  exchange  the  positions  of  the  source 
or  receiver.  But  the  monostatic  reverberation  intensities,  obtained  at  two  ter¬ 
minations  with  a  depth  of  Hi  or  H2,  would  not  be  reciprocal. 


^fli(’')  _  1^1  )N 


(17) 


where  N  is  the  angular  index  of  bottom  backscattering  defined  as 

BS{9)  =  10  logiQ  {y.Sin^B)  (18) 


If  Hi  =  50m,  H2  —  200m,  iV  =  2,  then  the  reverberation  interference  at 
shallower  area  would  be  VldB  less. 
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Introduction 

Strong  bottom  interaction,  variable  media  and  multi- 
path  propagation  make  shallow  water  acoustics  very  chal- 
lenging.  One  of  unexplained  phenomenon  for  shallow  wa- 
ter  sound  propagation,  in  the  contrast  with  the  so-called  *'  " 
the  optimum  propagation  frequency,  is  that  high  acoustic  j  *  * 
propagation  loss  over  some  frequency  range  is  frequently 
reported^~^.  Weston  attributes  such  a  phenomenon  to  fish 
activities^ .  In  some  case  it  is  associated  with  sediment  shear- 
wave  resonances  within  the  layer  of  sediment^.  Fig.  1  shows  v 
a  frequency  response  for  shallogv-water  sound  propagation  ^ 
in  the  summer  (with  a  strong  thermocline  shown  in  Fig.  * 

2)  in  the  Yellow  Sea  off  China,  obtained  by  Zhou  and  his  ' 
group  at  the  lAAS^.  Around  600  Hz  and  above  1100  Hz  ^ 
the  transmission  loss  is  abnormally  large.  Several  years  of  ' 
observations,  all  made  in  August  at  the  same  area  with  a  ! 

j  ^ 

similar  strong  thermocline,  have  shown  that  the  frequency  ;  ■>.  ^ 
responses  of  sound  propagation  is  often  a  strong  function  ' 
of  both  time  and  propagation  direction.  The  variability  of  ^  m 
the  high  loss  frequency  range  and  the  apparent  anisotropic  .'  ^ 
characteristics  of  sound  propagation  occur  with  no  obvious  ,  * '  ^ 
explanation  in  terms  of  fish  swim  bladder  resonance.  For  the  |  w>. 
strong  thermocline  shown  in  Fig.  2,  if  both  sound  source  and  I  7 
receiver  are  located  beneath  the  thermocline,  the  sea  surface 
influence  on  the  long-range  sound  field(or  lower  modes)  is 
negligible.  While  it  might  be  possible  that  one  could  use  the 
I  Weston  fish  model(or  some  other  model)  to  quantitatively 
match  the  experimental  data,  we  consider  the  possibility 
I  that  the  acoustic  normal-mode  conversion  caused  by  inter¬ 
nal  soliton  or  bottom  relief  packets  may  also  explain  the 
data. 

Ocean  models 

i.  Internal  Moliton  packets.  Our  numerical  results 
have  shown^  that  the  acoustic  mode-coupling  induced  by  | 
internal  wave  packets  with  a  gated  sine  waveform  (model  I)  ; 
exhibits  frequency,  soliton  wavelength  and  packet  length  res-  | 
onances.  The  interaction  between  the  acoustic  waves  and  in¬ 
ternal  wave  packets  could  provide  a  plausible  explanation  for  j 
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i  the  observed  anomalous  propagation  in  the  summer.  Now  . 

I  we  extend  the  model  to  consider  the  influence  of  the  cbarac-  | 
teristic  properties  of  more  realistic  oceanic  soliton  packets. 

Solitons  in  shallow  water  (/i/Ai  <<  1)  are  described 
to  first  order  in  wave  amplitude  by  the  Korteweg-de  Vries 
(KdV)  equation:  At  -h  QtAi  +  pAA,  +  6Axxx  =  0  .  A  so-  | 
lution  to  this  equation  is  A(x  —  ct)  »  o  *  sec/i*(S^).  The  i 
density  (temperature)  profiles  corresponded  to  Fig.  2  is  very  | 
reasonably  described  by  a  two-layer  fluid  model.  The  soli-  | 
ton  phase  speed  c  and  wavelength  Ai(and  scale  parameter  j 
L)  are  then  easy  to  obtain.  When  the  initial  solitary  wave  j 
propagates,  it  often  evolves  more  rank-ordered  solitons  and  | 
exhibits  clear  nonlinear  features.  Instead  of  the  previous 
simple  gated  sine  function,  three  packets  with  typical  diar- 
acteristics  of  internal  solitons  are  used;  the  classical  sech^ 
soliton  shape  w'hich  decreases  in  wavelength  and  amplitude  | 
toward  the  rear  of  the  packet(Model  H,  shown  in  Fig.  2  ). 

£.  Bottom  relief  group.  For  simplicity,  we  assume 
that  an  undulating  seabed  can  be  expressed  by  a  half  sine 
function,  shown  in  Fig.  2(Model  HI). 


I  The  presence  of  internal  soliton  or  bottom  relief  packets 
I  makes  the  environmental  parameters  range  dependent.  The 
,  parabolic  equation  (PE)  model  is  used  to  numerically  sim- 
{  ulate  the  effect  of  internal  soliton  packets  or  bottom  relief 
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groups  on  sound  propagation  in  shallow  water.  We  arbitrar¬ 
ily  put  three  packets  of  solitons  or  bottom  relief  located  at 

5  km,  15  km  and  25  km  along  propagation  track.  The  PE 
metbod^~^  is  used  to  calculate  the  frequency  response  of 
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acoustic  transmission  loss.  When  analyzing  the  characteris- 
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tics  of  acoustic  mode-coupling,  we  put  just  a  single  group 
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.  of  soUtons  or  bottom  relief  at  15  km  and  use  the  first  nor- 
,  mal  mode  alone  as  the  initial  input  field  to  the  PE  code 
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(IFD).  The  Acoustic  field  obtoined  using  P£  model  is  de¬ 
composed  into  the  normal  modes.  The  results  are  compared 
with  mode-coupling  theory^~^  in  order  to  improve  our  un¬ 
derstanding  of  how  mode  conversion  acts  to  become  an  im¬ 
portant  mechanism  of  acoustic  attenuation  in  shallow  water. 


1.  Similar  to  model  I  and  model  HI,  internal  virave  . 
padcets  with  the  classic  characteristic  properties  of  soliton  | 
could  also  cause  abnormally  large  attenuation  for  acoustic  ' 
propagation  in  shallow  water.  The  differences  between  the 
results  with  and  without  the  soliton  packets  is  shown  in  Fig. 

1  by  a  a  a  . 

2.  The  interaction  of  sound  waves  with  soliton  or  bot¬ 
tom  relief  group  exhibits  (signal)  frequency,  (radial  inho¬ 
mogeneity)  wavelength  and  packet  length  resonances.  As  ; 
an  example,  the  bottom  relief  wavelength  and  group  length 
resonances  are  shown  in  Fig.  3  (/  ■>  630.H’s,  r  m  30km).  i 

3.  For  a  shallow  water  with  a  strong  thermocline  the 
apparent  resonance  interaction  of  sound  waves  with  three 
ocean  models  could  only  occur  under  a  specific  circumstance: 
when  the  acoustic  mode  coupling  caused  by  environmental 
parameter  variation  transfers  a  significant  amount  energy 
from  lower  mode  into  higher-order  modes  with  much  larger 
attenuation  rate.  For  example,  the  PE  field  at  18km  is  de¬ 
composed  to  normal  modes  with  the  results  shown  in  Fig. 
4(f  ^  630 ffz,r  m  18km).  Apparently,  after  the  interaction 
of  the  first  mode  with  a  packet  at  15  km  a  significant  amount 
energy  has  been  coupled  into  higher-order  modes. 

4.  The  ’fresonacelike”  behavior  of  transmission  loss  pre¬ 
dicted  by  the  PE  analysis  is  consistent  with  noode  coupling 
theory.  Significant  energy  transfer  will  occur  between  mode 
m  and  n  if  kinj,  «  k,„  -  ik,,  here  k^^  is  the  wave  number  of 
the  spectrum  peak  of  the  radial  inhomogeneity.  The  mode 
coupling  is  a  periodic  (resonance)  function  of  the  soliton  (or  j 
relief)  group  length. 

5.  The  main  intention  of  this  work  is  to  show  what 
would  happen  to  sound  propagation  if  soliton  or  bottom  re-  ' 
lief  packets  are  present,  not  to  solve  an  inversion  problem 
for  internal  waves  or  to  determine  sound  propagation  for 
a  specific  area.  It  is  expected  that  any  abnormal  attenua-  ' 
tion  for  sound  propagation  in  a  real  shallow  wato'  environ-  ; 
ment  might  be  caused  by  ai^  combination  of  several  differ- 
ent  mechanisms,  internal  wave(or  bottom  relief)  interaction 
must  be  considered  as  one  such  mechanism. 
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Naturally  occurring  internal  solitary  wave  trains  (solitons)  have  often  been  observed  in  the 
coastal  zone,  but  no  reported  measurements  of  such  solitary  waves  include  low-frequency  long- 
range  sound  propagation  data.  In  this  paper,  the  possibility  that  internal  waves  are  responsible 
for  the  anomalous  frequency  response  of  shallow-water  sound  propagation  observed  in  the 
summer  is  investigated.  The  observed  transmission  loss  is  strongly  time  dependent,  anisotropic 
and  sometimes  exhibits  an  abnormally  large  attenuation  over  some  frequency  range.  The 
parabolic  equation  (PE)  model  is  used  to  numerically  simulate  the  effect  of  internal  wave 
packets  on  low-frequency  sound  propagation  in  shallow  water  when  there  is  a  strong 
thermocline.  It  is  found  that  acoustic  transmission  loss  is  sensitive  to  the  signal  frequency  and 
is  a  “resonancelike”  function  of  the  soliton  wavelength  and  packet  length.  The  strong 
interaction  between  acoustic  waves  and  internal  waves,  together  with  the  known  characteristics 
of  internal  waves  in  the  coastal  zone,  provides  a  plausible  explanation  for  the  observed 
anomalous  sound  propagation  in  the  summer.  By  decomposing  the  acoustic  field  obtained 
from  the  PE  code  into  normal  modes,  it  is  shown  that  the  abnormally  large  transmission 
attenuation  is  caused  by  “acoustic  mode-coupling”  loss  due  to  the  interaction  with  the  internal 
waves.  It  is  also  shown  that  the  “resonancelike”  behavior  of  transmission  loss  predicted  by  the 
PE  analysis  is  consistent  with  mode  coupling  theory.  As  an  inverse  problem,  low-frequency 
acoustic  measurements  could  be  a  potential  tool  for  remote-sensing  of  internal  wave  activity  in 
the  coastal  zone. 

PACS  numbers:  43.30.Bp.  43.30.Ft,  43.23.Cb 


INTROOUCTiON 

Naturally  occurring  internal  wave  packets  have  often 
been  observed  in  the  coastal  zone,  especially  in  the  summer. 
The  mechanism  for  the  generation  of  these  nonlinear  inter¬ 
nal  waves  have  been  widely  investigated  in  the  geophysics 
and  fluid  mechanics  community.  Unfortunately,  however, 
no  reported  measurements  of  such  solitary  waves  include 
low-frequency  long-range  sound  propagation  data.  The 
acoustic  community  has  paid  little  attention  to  the  effect  of 
solitary  waves  on  sound  propagation  with  the  exception  of 
the  work  of  Baxter  and  Orr  that  was  based  on  ray  theory  and 
calculated  the  influence  of  an  oceanic  internal  wave  packet 
on  short-range  (high-frequency)  sound  propagation.'  Ex¬ 
periments,  conducted  by  Zhou  and  his  group  at  Institute  of 
Acoustics  of  the  Chinese  Academy  of  Sciences  in  Beijing 
over  a  four-year  period  at  the  same  area  of  Yellow  Sea,  have 
shown  that  the  frequency  response  of  shallow  water  sound 
propagation  in  the  summer  is  a  strong  function  of  time  and 
propagation  direction,  and  sometimes  exhibits  an  abnormal¬ 
ly  large  attenuation  over  some  frequency  range.  A  part  of 
these  results  was  reported  before, but  it  cannot  be  ex¬ 
plained  by  a  conventional  range-independent  model  of 
sound  propagation  using  reasonable  bottom  acoustic  param¬ 
eters  and  an  average  sound-speed  profile. 

In  this  paper,  we  investigate  the  possibility  that  the 
anomalous  propagation  results  -tre  due  to  the  presence  of 


internal  waves.  First,  we  briefly  discuss  the  characteristics  of 
internal  wave  packets  in  the  coastal  zone.  We  then  review  the 
aforementioned  experimental  results  of  Zhou  et  al.  that  ex¬ 
hibited  the  anomalous  frequency  response.  In  Sec.  Ill,  we 
hypothesize  that  this  anomalous,  anisotropic  frequency  re¬ 
sponse  is  caused  by  the  influence  of  internal  wave  packets. 
We  support  the  hypothesis  with  numerical  simulation  re¬ 
sults  obtained  using  the  parabolic  equation  (PE)  propaga¬ 
tion  model.  In  Sec.  IV,  we  decompose  the  acoustic  field  ob¬ 
tained  by  using  the  PE  model  into  the  normal  modes,  and 
show  that  the  abnormally  large  transmission  loss  which  oc¬ 
curs  over  certain  frequency  ranges  is  due  to  “mode-cou¬ 
pling”  loss  induced  by  the  internal  wave  packets.  In  Sec.  V, 
we  show  that  “resonancelike”  behavior  of  the  attenuation  is 
consistent  with  mode  coupling  theory. 

I.  CHARACTERISTICS  OF  INTERNAL  WAVE  PACKETS 
IN  THE  COASTAL  ZONE 

Internal  waves  have  been  observed  almost  ever;^ere 
in  the  ocean.*  In  the  open  ocean  they  are  best  described  as  a 
stochastic  phenomenon  with  a  broadband  frequency  wave 
number  spectrum.’  However,  analysis  of  extensive  data  on 
internal  waves  in  the  coastal  zone,**’’  has  shown  that  these 
waves  exhibit  the  properties  of  solitons.  The  experimental 
data  includes:  current  and  temperature  measurements:  ver¬ 
tical  profiles  from  (ITD.  XBT,  and  acoustic  echo  sounding 
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devices;  ship’s  radar  and  satellite  (or  space  shuttle)  images 
obtained  at  optical  and  radar  frequencies.  We  are  primarily 
interested  in  the  effect  of  internal  solitons  on  long-range 
sound  propagation  in  the  coastal  zone,  and  not  the  mecha¬ 
nisms  for  the  generation  or  propagation  of  internal  wave 
packets  ( which  has  been  widely  investigated  in  the  geophys¬ 
ics  and  fluid  mechanics  communities).  We  thus  limit  our 
discussion  of  internal  waves  to  the  following  summary  of  the 
relevant  characteristics  of  internal  wave  packets  in  the  coas¬ 
tal  zone.  ■  ■ 

( 1 )  Internal  waves  in  shallow  water  are  frequently  ob¬ 
served  in  deterministic  groups  (wave  packets)  with  well- 
defined  wavelengths  that  are  describable  as  sohtary  waves 
(or  solitons).  These  waves  are  usually  observed  in  summer 
when  they  are  trapped  in  a  strong  and  shallow  seasonal  ther- 
mocline. 

(2)  Solitons  in  shallow  water  (A  //I,  <  1 )  are  described 
to  first  order  in  wave  amplitude  by  the  Korteweg-de  Vries 
(KdV)  equation:  A,  +  CgA,  -h/iAA,  +  Syi^  =0.  The 
solitary  waveforms  often  have  a  “sech^”  profile.  When  an 
initial  solitary  wave  propagates,  it  often  evolves  more  rank- 
ordered  solitons  and  exhibits  clear  nonlinear  and  dispersive 
features  such  as  higher-than-linear  group  velocity,  and  a  de¬ 
crease  in  wavelength  and  amplitude  toward  the  rear  of  the 
packet. 

( 3 )  The  number  of  wave  packets  are  highly  correlated 
with  the  strength  of  the  local  tides:  The  maximum  number 
occurs  during  spring  tides  and  the  minimum  number  occurs 
during  neap  tides.  Whether  or  not  they  occur  is  critically 
dependent  on  the  structure  of  the  background  shear  and 
stratification  profiles. 

( 4 )  Solitary  wave  packets  propagate  shoreward,  and  are 
often  generated  by  a  tidally  driven  flow  over  sills,  continental 
shelf  edges,  or  other  major  variation  in  underwater  topo¬ 
graphy. 

( 5 )  Because  of  the  shallowness  of  the  summer  seasonal 
thermocline  and  the  large  amplitude  of  the  coastal  internal 
waves,  strong  surface  expressions  of  solitons  have  been  ob¬ 
served  with  a  variety  of  remote  sensors,  including  photo¬ 
graphs  and  synthetic  apenure  radar  (SAR)  from  satellites 
and  space  shuttles. 

Naturally  occurring  large-amplitude  internal  solitons 


have  been  reported  in  many  coastal  zones  of  the  world  such 
as:  the  Massachusetts  Bay,*"*  "  the  New  York  Bight, 
Gulf  of  California.'®  Andaman  Sea  offshore  Thailand,’" 
the  Australian  North  West  Shelf,”  the  Sulu  Sea  between  the 
Philippines  and  Borneo,”  ”  off  the  coast  of  Portugal,'^  ” 
off  Hiunan  Island  in  the  South  China  Sea  and  off  the  Strait  of 
Gibraltar  in  Alboran  Sea,^’  the  Scotian  Shelf  off  Nova  Sco¬ 
tia.”  the  Celtic  Sea.”  and  so  on.  They  have  also  been  ob¬ 
served  in  lakes.  As  an  example,  an  excellent  record  of  an 
internal  wave  packet  is  shown  in  Fig.  1.  It  was  obtained  by 
Orr  using  high-frequency  acoustic  scattering  in  the  Massa¬ 
chusetts  Bay.’  *  The  period  of  the  waves  is  Doppler  shifted 
by  the  ship’s  speed  of  about  2.3  kn.  The  seasonal  thermocline 
is  displaced  by  30  m  (arrow  I )  and  the  stratified  point  scat- 
terers  ( zooplankton? )  at  30  to  40  m  are  displaced  more  than 
20  m  ( arrows  2  and  3 ) .  The  heavy  acoustic  backscattering  in 
the  vicinity  of  arrow  4  and  extending  in  an  oscillatory  pat¬ 
tern  throughout  the  figure  is  possibly  caused  by  turbulent 
mixing  in  the  thermocline.  The  heavy  scattering  near  73  m  is 
probably  caused  by  euphausiid  and  mysid  shrimp.  Figure  2. 
taken  from  the  work  of  (Dsbome  and  Burch,’  shows  how 
rough  and  smooth  bands  on  the  water  surface  can  be  caused 
by  internal  solitons.  Such  surface  expressions  have  been  de¬ 
tected  by  ship,  satellite,  and  space  shuttle  in  many  areas.  For 
example.  Fig.  3,  taken  from  the  work  of  Liu,’ '  is  a  line  draw¬ 
ing  of  the  internal  wave  packeu  in  the  New  York  Bight, 
observed  from  Seasat  satellite,  which  clearly  illustrates  the 
anisotropic  characteristics  of  internal  waves  in  the  coasul 
zone.  There  are  numerous  wave  groups,  all  propagating 
shoreward.  Each  group  consists  of  many  solitons.  llie  num¬ 
ber  of  groups  often  depends  on  the  local  tide  strength  and  the 
density  profile  of  water  column,  i.e.,  it  is  a  function  of  time. 

The  characteristics  of  internal  wave  packets  in  the  coas- 
*  tal  zone  summarized  here,  will  be  helpful  in  explaining  the 
anomalous  experimental  results,  introduced  in  the  next  sec¬ 
tion. 

II.  FREQUENCY  RESPONSE  OF  SHALLOW-WATER 
SOUND  PROPAGATION  IN  THE  SUMMER 

Zhou  and  his  group  at  Institute  of  Acoustics.  Academia 
Sinica  (lAAS)  have  measured  the  frequency  response  of 
shallow-water  sound  propagation  under  the  condition  of  a 


FIG.  I.  A  200.kHz  acoustic  record  of  an  in- 
Mmai  wave  packet  otaerved  in  the  Massa- 
chuaetts  Bay. 
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FIG.  2.  A  ikclch  of  an  inicnul  Mlilon  in  ■  t«o-layer  fluid  model  showing 
the  fluid-paniek  vckxnia  tnd  me  surface  eapreaaioo. 

Strong  theimocline  in  the  Yellow  Sea  off  China.  A  great  deal 
of  dau  was  collected  over  several  years  in  the  same  area, 
(which  had  a  charaaeristically  flat  seabottom  with  high¬ 
speed  sediments).  The  data  included  measurements  of 
sound  propagation,  long-range  reverberation,  sound  field 
spatial  coherence  and  utilized  normal-mode  spatial  filtering 
techniques.  The  acoustic  parameters  of  the  bottom  for  this 
area  were  obtained  by  wideband  acoustic  measure¬ 
ments.’*'^^  Generally  speaking,  the  experimental  data  fiu 
theoretical  predictions  very  well.  For  example,  transmission 
losses  as  a  Unction  of  frequency  for  winter,  late  spring,  and 
sometimes  for  summer  were  correctly  predicted.  However,  it 
was  found  that,  in  the  summer,  even  along  the  same  experi- 
menul  track,  with  similar  averaged  sound-speed  profiles  the 
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FIG  3  Ltac  drawmg  indieatmg  mtcrnal  wivt  packMs  m  the  New  Yorii 
Bigbi,  otncnwi  firom  the  S«awi  iMclUit. 
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frequency  response  of  sound  transmission  could  often  be 
very  different  for  different  years. 

Figure  4  shows  power  spectra  that  are  obtained  by  aver¬ 
aging  several  explosive  propagation  signals.  The  upper  curve 
is  the  power  spectrum  at  a  propagation  range  of  0.5  km.  the 
lower  one  is  at  a  range  of  28  km.  The  difference  between  the 
two  curves  is  a  measure  of  the  sound  transmission  loss  be¬ 
tween  two  receiving  points,  i.e.,  the  frequency  response  of 
sound  propagation,  ^tween  KX)  and  1 100  Hz,  especially 
around  600  Hz,  the  transmission  loss  is  abnormally  large. 
The  measured  sound-speed  profiles  at  the  receiving  ship  are 
shown  in  Fig.  S.  Both  source  and  receiver  were  located  below 
the  theimocline.  The  anomalous  transmission  loss  cannot  be 
explained  by  conventional  models  of  sound  propagation 
(with  a  fill]  complement  of  bottom  acoustic  parameters  and 
an  averaged  sound-speed  profile  of  water  column ) .  In  order 
to  explain  this  phenomenon,  2niou  and  his  colleagues  at 
lAAS  made  observations,  every  August  for  4  years  at  same 
area.  The  abnormal  attenuation  frequency  range  varied  with 
time.  In  Fig.  4,  it  occurs  at  around  600  Hz,  but  for  different 
years  or  different  propagation  directions,  it  has  been  ob¬ 
served  around  500, 1 2()0,  or  1 600  Hz.  Sometimes,  there  is  no 
apparent  abnormal  attenuation  frequency  range. 

In  another  experiment,  the  frequency  response  at  a  fixed 
range  for  six  different  radial  directions  was  obtained.  The 
resultt  are  shown  in  Fig.  6.  The  transmission  distance  was 
kept  constant  at  28  km,  by  moving  the  source  ship  along  one 
quadrant  of  a  circle  centered  at  the  receiving  ship  Each 
curve  represents  an  averaged  value  obtained  from  several 
explosive  signals.  Transmission  loss  is  obviously  a  strong 
function  of  propagation  direction.  For  different  propagation 
directions,  at  some  frequencies,  the  sound  intensity  varied  as 
much  as  25  dB!  However,  in  another  experiment,  at  the  same 


FREQUENCY  (Hz) 

FIO.  4.  EspIciiveiisiialpowtfiptctnMn  far  ihaltow-wawr  sound  prop«ei- 
Bon  in  the  Mnmntr  which  exUbiu  m  abnermally  Ivtc  aitcnuuioo  around 

400 Hz  ( axpemnanui  dau ) .  Both  iourcc  and  racavCT  ate  locaiad  below  I  he 

tharmochna. Cawpuutional ramlu (see Sac. Ill  B):(l)whhautpacket-- 
(2)  with  three packen  O O O. 

Zheu  era/  Sound  and  auemaisoMons  2044 


SOUND  VELOOTY  (m^) 
uoo  isio  isao  ifiao  1640 
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place,  in  a  different  year,  two  orthogonal  propagation  direc¬ 
tions  were  observed  to  have  almost  identical  frequency  re¬ 
sponses.  It  was  not  at  all  clear  why.  when  there  is  a  strong 
Ihermocline,  the  frequency  response  of  shallow-water  sound 
propagation  in  the  experimental  area  was  often  observed  to 
be  a  strong  function  of  both  time  and  propagation  direction, 
and  sometimes  exhibited  abnormally  large  attenuation  over 
some  frequency  range. 

Using  ray  theory  and  expenmental  contours  of  sound 
speed  versus  range  and  depth.  Baxter  and  Orr'  calculated 
the  influence  of  an  internal  wave  packet  on  short  range 
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FIG-  7.  A  record  of  temperature  venus  tune  over  a  penod  of  42  min  at  a 
ftaed  depth  (near  the  thermocline). 


(high-frequency)  sound  propagation.  The  calculations 
showed  that  at  a  fixed  position  the  intensity  of  sound  field 
could  vary  by  as  much  as  20  dB,  relative  to  a  case  without  the 
presence  of  the  internal  wave  packet.  We  hypothesized  that 
internal  wave  packets  might  have  a  strong  influence  on  low- 
frequency  long-range  sound  propagation  in  shallow  water  as 
well  and  could  possibly  explain  Zhou's  Yellow  Sea  propaga¬ 
tion  dau.  In  the  next  section,  we  support  this  hypothesis 
with  numerical  simulation  results  obtained  using  the  para¬ 
bolic  equation  method. 

IIL  NUMERICAL  SIMULATION  OF  THE  INFLUENCE  OF 
INTERNAL  WAVE  PACKETS  ON  SOUND 
PROPAGATION 

A.  OcMiiie  modal 

The  presence  of  internal  waves  makes  the  sound-speed 
profile  of  the  water  range  dependent.  Unfortunately,  in 
Zhou's  experiments  there  were  no  accompanying  systematic 
measurements  of  the  internal  wave  field;  and  only  the  sound- 
speed  profile  at  the  receiver  was  measured.  During  the  time 
^  period  over  which  the  frequency  response  shown  in  Fig.  4 
was  obtained,  it  was  found  that  the  thermocline  depth  at 
receiving  ship  location  did  change  with  time.  Figure  5  shows 
three  sound-speed  profiles  that  were  obtained  at  different 
times. 

During  the  experiment.  Zhang'*  measured  the  tem¬ 
perature  fluctuation  as  a  function  of  time  at  the  receiving 
ship.  A  record  of  temperature  fluctuations  over  a  period  of 
42  min  at  a  fixed  depth  ( around  the  thermocline )  is  shown  in 
IRg.  7.  The  temperature  fluctuated  between  about  13  *C  and 
27  *C.  The  possible  presence  of  several  individual  solitary 
waves  is  indicated  by  temperature  peaks  due  to  the  depres¬ 
sion  of  the  thermocline.  Figures  S  and  7  give  at  least  some 
evidence  that  there  was  solitary  wave  activity  in  the  area  of 
acoustic  experiments. 

Due  to  the  fact  that  we  have  no  data  concerning  specific 
characteristics  of  internal  waves  at  the  experimental  site,  for 
simplicity,  following  Lee’s  three-layer  model  of  internal 
wave,**  we  assume  that  the  internal  wave  packet  can  be  ex¬ 
pressed  by  a  gated  sine  function  as  shown  in  Fig.  8(b).  In 
this  figure,  if  a  packet  begins  at  /-g,  then,  for  r,  -f- 

Z,  «  14.0  —  2.0  m(2irr/A, ),  ( 1 ) 

Z,  *  19.0— 5.0 5in(2iiT/A,).  (2) 

The  idealized  sound-speed  profile  in  the  absence  of  in¬ 
ternal  waves  for  the  numerical  simulation  is  shown  in  Fig. 
8(a).  We  call  A,  the  solium  wavelength,  and  the  packet 
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FIG.  8.  Simplified  oceanic  model  for  an  inlemal  wave  packet. 


length.  Referring  back  to  Fig.  1,  we  see  that  an  almost  sinu¬ 
soidal  interface  exists  between  the  wanner  surface  layer  and 
the  underlying  cooler  water.  The  waveform  of  the  soliton  in 
Fig.  5  of  Ref.  1  is  also  almost  sinusoidal.  The  simplified, 
sinusoidal  model,  is  thus  seen  to  be  not  physically  unreason¬ 
able.  The  main  purpose  of  this  work  is  to  provide  the  first 
physical  insight  into  the  possible  effects  of  internal  wave 
packets  on  wideband  low-frequency  sound  propagation  in 
shallow  water.  It  is  not  to  predict  the  characteristics  of  inter¬ 
nal  waves  for  a  specific  area.  Our  simplification  should  not 
alter  the  qualitative  results  which  follow. 

The  averaged  values  for  bottom  sound  speed  (c» )  and 
bottom  acoustic  attenuation  (a^ )  for  the  experimental  area 
are  given  by^*-“ 

c,/c„.  =  1.056.  .(3) 

a,  =  0.34/' dB/m,  (4) 

where /is  the  frequency  in  kHz. 

B.  Comparison  of  numerical  frequency  response  wl'.h 
experimental  data 

Numerical  simulation  results  obtained  using  the  para¬ 
bolic  equation  (PE)  method  (IFD^'"*’  code),  show  that 
internal  wave  packets  can  significantly  change  low-frequen¬ 
cy  transmission  loss.  The  frequency  response  and,  in  particu¬ 
lar,  the  abnormal  attenuation  frequency  range  are  sensitive 
functions  of  the  parameters  of  the  soliton  wave  packets. 

For  example,  using  the  simplified  ocean  model  given  in 
Fig.  8  and  Eqs.  ( 1  )-(4),  we  put  three  internal  wave  packets 
located  at  5,  IS,  and  25  km,  along  propagation  track  with 
each  packet  consisting  of  six  solitons  each  with  a  wavelength 
Ai  of  235  m  (£,  =  =  1410  m).  The  transmission  loss  as 

a  function  of  range  for  four  frequencies  (with  a  receiver 
depth  of  32  m  and  a  source  depth  of  25  m )  is  shown  in  Fig.  9. 
The  difference  between  the  results  with  and  without  the  soli¬ 
ton  packett  present  is  small  at  100  or  lOOO  Hz,  but  at  6(X)  Hz 
the  difference  at  30  km  is  as  much  as  25  dB.  The  most  inter¬ 
esting  and  encouraging  result  is  that,  in  this  case  the  numeri¬ 
cal  transmission  loss  as  a  function  of  frequency  (shown  in 
Fig.  4  by  the  circles)  fits  the  experimental  results  quite  well. 
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As  a  consequence  of  the  characteristics  of  internal  wave 
packets  described  in  Sec.  I  (propagation  in  groups,  well- 
defined  wavelengths,  high  correlation  with  tides,  shoreward 
propagation),  the  projection  of  the  solitons  in  different 
acoustic  propagation  directions  and  for  different  times  will 
be  different.  If  we  change  the  soliton  wavelength  in  the  three 
packets  to  235,  3(X},  350, 400,  and  500  m,  we  would  get  fre¬ 
quency  responses  as  shown  in  Fig.  10,  that  are  similar  to  the 
experimental  results  for  different  propagation  directions 
shown  in  Fig.  6.  Hence,  interaction  between  acoustic  waves 
and  internal  wave  packets  is  consistent  with  Zhou's  experi- 
menul  results  for  shallow-water  sound  propagation  in  the 
summer. 

Why  does  the  abnormally  large  attenuation  in  Fig.  4 
occur  around  630  Hz  for  internal  wave  packets  consisting  of 
six  solitons  with  a  wavelength  of  235  m?  Why  do  internal 
wave  packets  with  the  same  parameters  have  a  much  smaller 
effect  on  sound  propagation  at  frequencies  around  300  or 
1000  Hz?  In  the  next  section,  we  will  show  that  these  are  due 
to  “resonances”  in  the  acoustic  mode-coupling  induced  by 
the  internal  wave  packets. 

IV.  RESONANCE  EFFECTS 
A.  Frequuncy  resonance 

For  our  numerical  model  of  shallow  water,  with  a  water 
depth  of  38  m,  mode  stripping  caused  by  seabottom  attenu¬ 
ation,  results  in  only  a  few  modes  still  being  present  at  a 
distance  of  15  km.  If  the  source  and  receiver  are  located 
below  the  thermocline,  for  a  realistic  sea  bottom,  the  first 
mode  dominates  the  sound  field  at  long-range  over  the  fre¬ 
quency  range  of  interest.  In  order  to  isolate  the  effect  of  in¬ 
ternal  wave  packets  on  the  mode  coupling  among  acoustic 
normal  modes,  we  put  a  single  packet  conristing  of  six  soli- 
tons  with  a  wavelength  of  235  m  at  15  km  and  use  the  first 
normal  mode  akme  as  the  initial  input  field  to  the  parabolic 
equation  code.  Thus,  prior  to  interaction  with  the  packet,  we 
assume  only  one  mode  (the  first)  is  present;  higher-order 
modes  are  generated  as  the  first  mode  propagates  through 
the  internal  wave  packet.  The  depth  distribution  function  of 
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FIG.  9.  Comparison  between  transmission  lo»  without  internal  waves  (•--)  and  with  three  packets  located  at  S.  IS.  and  2S  km  ( — ).  L,  «  1410  m. 
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sound  pressure  amplitude  at  1 3, 20,  and  30  km  obtained  by  ever,  for  630  Hz,  the  depth  distribution  shape  is  quite  differ- 

PE  method  are  shown  in  Fig.  11.  At  a  range  of  13  km  (before  ent  from  the  first  mode,  due  to  resonance  scattering,  and  the 

the  wave  has  interacted  with  the  internal  wave  packet ),  the  amplitude  is  much  smaller. 

shape  ofthe  first  mode  is  maintained  at  all  fiequencies.  After  The  mode  coupling  due  to  interaction  with  internal 
interaction  with  the  internal  wave  packet,  at  300  and  1000  wave  packets  can  be  easily  calculated.  Neglecting  the 

Hz,  the  energy  is  still  predominantly  in  the  first  mode.  How-  evanescent  modes,  which  decay  at  large  range,  leaves  only  a 
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finite  number  of  propagating  modes  at  a  given  frequency. 
The  propagated  field  [d>pE  (r^)  ]  obtained  by  the  PE  meth¬ 
od  can  be  expanded  in  terms  of  a  set  of  local  modal  eigen- 
fimctions  [I/, (z)],  corresponding  to  the  sound  velocity 
profile  shown  in  Fig.  8(a),  as  follows: 

(5) 

m 

where  the  local  modal  eigenfunctions  U,  (z)  satisfy  the  ap¬ 
propriate  boundary  conditions  at  the  sea  surface  and  bottom 
and  the  differential  equation 

i^^[k\z)-k\]U,=Q,  (6) 

where  A;,  is  the  eigenvalue  ( modal  wave  number ) .  The  mod¬ 
al  eigenfunctions  also  satisfy  the  orthogonality  relation 

|p(z)f/,(z)£/„(z)rfz  =  5,.„.  (7) 

The  Kronecker  delta  on  the  right  side  is  unity  ifn  =  m  and  is 
zero  if  n^m. 

Multiplying  both  sides  of  Eq.  (5)  by  U„(2)  and  inte¬ 
grating  over  depth  with  help  of  Eq.  ( 7 )  ( mode  spatial  filter¬ 
ing),  we  obtain 
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FIG.  1 1 .  The  depth  distribution  function  of  sound-pressure  amplitude  for 

different  frequencies  at  a  range  of  IS  km  ( - ),  20  km  ( — ).  and  30  km 

(— )  obtained  using  the  PE  method.  PE  input  Held:  mode  1.  One  internal 
wave  packet  at  IS-I6.41  km  with  A,,  •  235  m. 
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FIG.  12.  The  retaiiveimpliiuda  of  different  modejK  1  *  km.  showing  frequency  resonance  of  mode  coupling.  PE  input  field:  mode  I.  One  packet  at  15-16.41 
km.  A,  w  235  m. 


Af^ir)  =  j  p(z)<Pr^{r^)U,(z)dz. 


After  the  interaction  of  the  first  mode  with  the  internal 
wave  packet,  at  selected  ranges  where  no  internal  wave  ex¬ 
ists,  the  relative  amplitude  of  higher-order  mode  to  the  first 
mode,  is  given  by 

^  _  Mpe.(>')I  _\Sp{z)^n(r^)U,{z)dz\ 

Mpei('’)|  \fp(z)<P„(r^)U,(z)dz\ 


^  /A,j,+A„y'^ 

A-  Ail) 

b 

A,i,^(f,Pi<PM(>’<2j)UJZj) 

\j- 1 

Jm  m  l  ' 

\m  I 

+  T  Pt^i{rj:j)U.{Zj)\ , 
y-M>i  / 


where  <!>«  ( r,2, )  and  )  are  the  real  part  and  the  imagi¬ 

nary  part  of  PE  field  in  the  yth  depth  increment,  Af  is  the 
number  of  layers  and  the  water-bottom  interface  occurs  be¬ 
tween  the  layer  m  and  layer  (m -(-  1).^‘  The  quantity,  is 
a  measure  of  the  energy  exchanged  between  the  first  mode 
and  higher-order  modes.  In  this  manner  using  a  normal¬ 
mode  computer  program,  a  PE  computer  program  and  the 
ocean  model  described  in  Sec.  Ill,  we  were  able  to  decom¬ 
pose  the  sound  field  obtained  from  PE  into  normal  modes. 
The  input  field  for  the  PE  code  is  the  first  normal  mode.  In 
the  absence  of  modal  coupling  only  this  mode  would  ever  be 
present.  Results  for  three  different  frequencies  are  shown  in 
Fig.  12.  The  plots  show  the  relative  modal  amplitudes  at  1 8 
km.  For  ranges  less  than  IS  km,  only  the  first  mode  is  pres¬ 
ent  As  indicated  in  Fig.  12,  beyond  IS  km  there  is  interac¬ 
tion  with  the  internal  wave  packet  and  at  630  Hz  a  signifi¬ 
cant  amount  ofenergy  has  been  transferred  from  mode  I  into 
higher-order  modes  that  will  attenuate  rapidly  thereafter. 
Relatively  little  modal  coupling  occurs  at  the  other  two  fre¬ 
quencies.  Table  I  shows  the  eigenvalue  Ic,  and  the  attenu¬ 
ation  rate  0^  (corrected  for  cylindrical  spreading)  for  the 
first  seven  normal  modes,  calculated  by  a  normal  mode  pro¬ 
gram  for  the  sound  velocity  profile  shown  in  Fig.  8(a)  and  a 
firequency  of  630  Hz. 
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TABLE  I.  The  ofcnvaliN  ,  tncniution  ntt  lad  afCDvmliw  diffmnce 
between  modes  I  end  •  for  the  tm  seven  noraul  modes, />•  630  Hz- 


Mode  It 

0,  (dBAm) 

k,  —  k. 

1 

2.633  594  3 

0.2434 

0 

2 

2.626  323  6 

0.9169 

0.009  2687 

3 

2.611  159  6 

1.7661 

0.024434  7 

4 

2.391  298  6 

2.4771 

0.044295  7 

5 

2.572  733  4 

1.3938 

0.062  8389 

6 

L363  417  2 

1.7853 

0.072  177  1 

7 

2.348  399  6 

2.3948 

0.087  194  7 

From  this  analysis,  it  is  evident  that  the  modal  coupling 
caused  by  internal  waves  can  sometimes  be  an  important  loss 
mechanism  for  sound  transmission  in  shallow  water  in  the 
summer.  At  300  and  1000  Hz,  the  mode-coupling  effect  is 
much  weaker.  Only  a  few  percent  of  the  wave  energy  is  cou¬ 
pled  into  high-order  modes.  At  630  Hz,  however,  a  signifi¬ 


cant  amount  ofenergy  has  been  transferred  from  mode  1  into 
higher-order  modes,  especially  mode  3  and  mode  4,  with 
much  larger  attenuation  rate.  The  mode-coupling  induced 
by  internal  wave  packets  exhibiu  a  frequency  resonance  ef¬ 
fect,  and  can  cause  abnormally  large  transmission  loss 
around  the  resonance  frequency.  The  abnormally  high  at¬ 
tenuation  observed  by  Zhou  et  ai  over  certain  frequency 
ranges  is  consistent  with  interaction  with  internal  waves.  In 
Sec.  V,  we  will  discuss  what  determines  the  “resonance  fre¬ 
quency.” 

B.  Sollton  wavolongtti  reaoruinee 

We  continue  to  consider  a  single  packet  located  at  a 
distance  of  IS  km.  We  now  fix  the  acoustic  frequency  at  630 
Hz  and  the  packet  length  (Lf)  at  1.4-l.S  km  and  calculate 
the  effect  of  soliton  wavelength  (Aj )  on  transmission  loss. 
Figure  13(a)  shows  the  transmission  loss  as  a  function  of 
soliton  wavelength  at  a  range  of  30  km  ( with  a  receiver  depth 


(•) 


(c) 


(b)  (d) 


FIG.  13.  (i)  Soliton  wavelength  reso¬ 
nance.  (b)-(d)  The  depth  distribution 
function  of  sound-pressure  amplitude  for 
different  soliton  wavelengths  at  15  km 
( - ),  20  km  ( — ).  and  30  km  (— )  ob¬ 

tained  by  PE  method.  PE  input  field: 
mode  1.  One  packet  at  15-16.41  km. 
/-630Hr. 
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FIO.  14.  The  relalive  amplitudes  of  different  modes  at  1 8  km,  showing  soliton  wavelength  resonance  of  mode  coupling.  PE  input  held:  mode  I  One  packet  at 
15-16.41  km./=630Hz. 


of  32  m  and  a  source  depth  of  23  m).  Transmission  loss 
exhibits  a  resonancelike  maximum  at  a  soliton  wavelength  of 
233  m.  For  other  soliton  wavelengths  between  120  and  330 
m,  the  difference  between  transmission  loss  with  and  with¬ 
out  the  presence  of  an  internal  wave  packet  is  not  significant. 
Taking  first  normal  mode  as  the  input  field  to  the  PE  code, 
we  obtained  the  depth  distribution  functions  shown  in  Fig. 
13(b)-(d).  After  interaction  with  the  packet,  for  a  soliton 
wavelength  of  233  m,  the  shape  is  very  different  than  that  of 
the  first  mode,  but  for  200  and  280  m,  the  shapes  are  very 
close  to  that  of  mode  1. 

At  a  distance  of  18  km.  we  again  decompose  the  PE 
calculated  sound  field  into  normal  modes  and  get  the  results 
for  three  different  soliton  wavelengths  shown  in  Fig.  14.  For 
a  soliton  wavelength  of  233  m.  a  significant  amount  of  ener¬ 
gy  has  been  transferred  from  mode  1  into  higher-order,  high¬ 
er-loss  modes.  However,  for  soliton  wavelengths  of  200  and 
280  m,  the  mode-coupling  effect  is  rather  wodt.  That  is  the 
mode  coupling  and,  hence,  the  loss  induced  by  internal  wave 
packets  exhibits  a  soliton  wavelength  resonance  effect. 

C.  Packat  lofiflth  raaonanca 

Once  again,  considering  a  single  packet  located  at  a 
range  of  13  km,  we  now  hold  acoustic  frequency  (630  Hz) 


and  soliton  wavelength  ( 233  m )  constant,  and  calculate  the 
effect  of  internal  wave  packet  length  (or,  equivalently,  the 
number  of  solitons)  on  transmission  loss.  Figure  13 (a) 
shows  the  transmission  loss  as  a  function  of  packet  length  at 
30  km  ( with  a  receiver  depth  of  32  m  and  a  source  depth  of 
23  m).  The  transmission  loss  is  a  periodic,  “resonancelike” 
function  of  the  packet  length  (Z,).  Taking  first  normal 
mode  as  the  input  field  to  the  PE  code,  we  obtain  the  depth 
distribution  functions  for  the  sound  field  for  three  packet 
lengths  shown  in  Fig.  13(b)-(d).  At  18  km,  after  interaction 
with  the  packet,  the  shape  of  the  depth  distribution  function 
fora  six-soliton  packet  is  very  different  from  that  of  the  first 
mode.  However,  for  a  longer  packet  consisting  of  1 1  solitons 
it  is  close  to  the  shape  of  the  first  mode  and  close  even  in 
amplitude  to  the  results  without  a  packet. 

At  a  distance  of  18  km,  we  again  decompose  the  PE 
sound  field  into  normal  modes,  and  get  the  results  for  two 
different  packet  lengths  shown  in  Fig.  16.  For  a  14l0-m 
packet  (six  23S-m  solitons),  a  significant  amount  of  energy 
has  been  transferred  from  mode  1  into  higher-order  modes. 
For  a  longer  packet  (2583  m,  i.e.,  for  1 1  solitons),  the  mode- 
coupling  effect  is  much  weaker.  That  is,  the  mode-coupling 
and  hence  the  loss  induced  by  internal  wave  packets  exhibits 
a  packet  length  resonance  effect. 
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FIG.  IS.  (a)  Packet  length  resonance,  (b)-(d)  The  depth  distribution 
function  of  sound-pressure  amplitude  for  different  packet  lengths  at  I S  km 

( - ),  20  km  (—),  and  30  km  (— )  obtained  by  PE  method.  PE  input  held: 

mode  I.  One  packet  beginning  from  IS  km./^i  630  Hz. 


V.  MODAL  COUPUNG  ANALYSIS 

Dozier  and  Tappert,^  and  McDaniel  and  McCam- 
mon^’  have  shown  that  the  exchange  of  energy  between 
acoustic  modes  n  and  m  induced  by  random  internal  waves 
in  the  deep  sea  or  lateral  seabed  inhomogeneities  in  shallow 
water  exhibits  a  wave  number  resonance  effect.  The  princi¬ 
pal  transfer  of  energy  will  occur  between  modes  whose  eigen¬ 
value  difference  equals  the  wave  number  of  the  spectrum 
peak  of  the  intern^  waves  or  sub-bottom  roughness.  This 
result  is  also  applicable  to  our  numerical  results  shown  in 
Figures  11-14.  Significant  energy  transfer  will  occur  be¬ 
tween  modes  m  and  n  if 

where  =  hr/Aj  is  the  wave  number  of  the  internal  wave. 
[  If  the  mode  couping  is  considered  to  be  a  quantum  mechan¬ 
ical  phonon-soliton  interaction,  Eq.  (12)  is  a  statement  of 
conservation  of  momentum.  The  corresponding  statement 
of  conservation  of  energy  a>,„,  =0,  —o)„  can  be  shown  to 
correspond,  macroscopically,  to  the  Doppler  shift  which  re¬ 
sults  from  scattering  from  the  moving  internal  wave.  ]  For 
our  numerical  model  of  internal  wave  packets  with  a  wave¬ 
length  of 235  m,  the  wavenumber  of  the  spectrum  peak  of  the 
internal  wave  k,„,  s;2tr/23S  »  0.026  737  0.  Ck>mpaTing  this 
value  of  A;,„,  with  the  eigenvalue  differences  between  mode  1 
and  mode  n  shown  in  Table  I,  it  is  found  that  k„,  is  almost 
equal  to  the  eigenvalue  difference  between  acoustic  normal 


FIO.  16.  The  relative  amplitudes  of  dilferent  modes  at  18  km.  showing  packet  length  resonance  of  mode  coupling.  PE  input  field;  mode  I.  One  packet 
begmiung  from  IS  km./w  630  Hz.  A,  •  23S  m. 
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1  and  3.  Thus  one  would  expect  a  significant  amount 
gy  to  be  transfered  from  mode  1  to  mode  3  ( but  not, 
mode2).  Moreover,  is  also  reasonably  close  to  the 

ices  between  the  eigenvalues  of  modes  3  and  4  and 
a  modes  4  and  6,  etc.  Hence,  one  would  expect  some 
X  of  energy  from  mode  1  to  mode  4  via  mode  3,  or 
node  1  to  mode  6  via  modes  3  and  4,  and  so  on.  These 
;  arguments  explain  why  a  significant  amount  of  ener- 
ransferred  from  mode  1  into  higher-order  modes,  espe- 
modes  3  and  4,  and  why  and  ^here  the  acoustic  wave 
sacy  and  soliton  wavelength  resonances  occur, 
lie  packet  length  resonance  eflTect  mentioned  in  Sec.  IV 
c  so  easily  explained,  and  it  is  still  not  clear  how  to 
ct  where  it  will  occur.  We  present  here  a  preliminary 
nation  for  the  phenomenon  given  by  S.  T.  McDaniel.^* 
E  backscattered  field  is  neglected,  the  equation  govem- 
lode  coupling  between  mth  and  nth  modes  due  to  range 
ndence  of  the  ocean  environment,  can  be  expressed  as 
Eq.  ( 14)  in  Ref.  37  for  noution] 

=  exp[i(k„ -k,)r].  (13) 


Let  a  resonant  coupling  condition  be 

=  -  2JCt  exp [  - 1( k„  -  k,  )r]  ( 14) 

I  Af„,  =  —  .  We  consider  only  two  modes,  modes  1 

1  m,  which  are  assumed  to  satisfy  the  resonant  coupling 
idition.  From  Eqs.  ( 13)  and  ( 14),  we  obtain 


-7-= 

dr 

(13) 

a 

II 

(16) 

,  ..2 

(17) 

If  all  of  the  energy  is  assumed  to  be  initially  in  the  first 
ide  (u,  =  1  and  =  0}  the  solution  to  Eq.  (17)  is 


u,  sscoslfiT^r),  (18) 

=sin(/rir).  (19) 

Equation  (19)  predicts  that  the  magnitude  of  the  excit- 
higher-order  modes  will  be  periodic  with  range.  The  reso- 
at  coupling  parameter  determines  for  what  packet 
igths  u„  has  a  maximum  value  (or  u,  has  a  minimum 
ue).  From  Eq.  ( 19),  the  packet  length  “resonance”  will 
;ur  when  =»(«-)•  l/2)ir.  From  the  first  packet 
igth  resonance  in  Fig.  lS(a),  we  have 
235x6a>4r/2  and  A^sLllAii;  from  the  second 
onance  length,  we  have  A^I,  »  x235x  17  *  3ir/2 

3  A4  at  1.18/km.  The  two  values  obtained  for  A^  areal- 
St  the  same  The  agreement  is  quite  satisfactory  consider- 
;  that  the  mode  coupling  occun^  between  more  than  just 
E>  modes  (modes  I  and  3)  and  A(.  would  be  changed  with 
!  coupling  strength.  The  mode  coupling  is  a  periodic  (res- 
ince)  (unction  of  the  internal  wave  packet  length,  and, 
Bce  acoustic  transmssion  loss  at  long-range  must  vary 
th  the  packet  length.  It  is  shown  that  the  “resonancelike" 


behavior  of  transmission  loss  predicted  by  the  PE  analysis  in 
this  paper  is  consistent  with  mode  coupling  theory. 


VI.  CONCLUSION 

In  summary,  ( 1)  we  have  briefly  reviewed  the  charac¬ 
teristics  of  naturally  occurring  internal  wave  packets  in  the 
coastal  zone  which  differ  in  many  respects  from  open  sea 
internal  waves.  We  demonstrate  that  such  internal  wave 
packets  can  have  a  strong  influence  on  low-frequency  long- 
range  sound  propagation  in  shallow  water.  (2)  The  mea¬ 
sured  frequency  response  of  sound  propagation  during  the 
summer  is  often  a  strong  function  of  both  time  and  propaga¬ 
tion  direction,  sometimes  the  sound  propagation  has  an  ab¬ 
normally  targe  attenuation  over  some  frequency  range.  (3) 
Numerical  calculations  have  shown  that  the  interaction  be¬ 
tween  the  acoustic  waves  and  internal  wave  packets  and  par¬ 
ticularly  resonance  effects  in  the  acoustic  mode-coupling  in¬ 
duced  by  internal  wave  packets  could  be  an  imporunt  loss 
mechanism  for  shallow-water  sound  propagation  in  the 
summer.  Modal  coupling  induced  by  internal  waves  could 
explain  the  anomalous  acoustic  experimental  results.  (4) 
The  fact  that  the  acoustic  mode-coupling  induced  by  inter¬ 
nal  wave  packets  exhibits  frequency,  soliton  wavelength,  and 
packet  length  resonances  [and  perhaps  existence  of  the 
Doppler  shift  alluded  to  in  Sec.  V,  as  well]  suggest  that  low- 
frequency  acoustic  mea''>rements  could  be  used  for  remote 
monitoring  of  internal  wave  activity  and  extracting  of  hy¬ 
drological  and  meteorological  characteristics  of  the  water 
mass  in  the  coastal  zone. 

The  analysis  and  numerical  calculations  presented  here 
are  based  on  a  simplified  oceanic  model.  Although  the  sim- 
,  plification  should  not  alter  the  qualitative  results  presented, 
a  more  refined  model  would  be  desirable  for  a  more  detailed 
comparison  with  observations.  It  is,  of  course,  also  possible 
that  some  seabottom  structure  (for  example,  a  surficial  layer 
of  low  speed),  mode-coupling  due  to  the  seabed  roughness 
or  sediment  inhomogeneities  and  fish  shoals  could  also  pro¬ 
duce  an  extra  acoustic  attenuation  over  some  frequency 
range.  It  would  thus  be  desirable  for  the  acoustic,  marine 
geology  and  geophysics  communities  including  remote- 
sensing  groups  to  work  together  to  conduct  joint  at-sea  ex¬ 
periments  at  a  specific  sea  area.  If  data  on  internal  wave 
activity  (solitary  or  random),  seabottom  parameters,  and 
sound  propagation  were  obtained  simultaneously  and  sys¬ 
tematically.  our  understanding  of  the  interaction  between 
internal  waves  and  sound  waves  in  the  coastal  zone  and  low 
frequency  acoustic  propagation  loss  in  shallow  water  would 
be  greatly  enhanced. 
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